Coherent Control of Quantum Chaotic Diffusion 



O 

o 
o 

(N 
O 

Q 

(N 



> 
O 
in 

(N 

o 
o 



9 L,: 



13 



Jiangbin Gong and Paul Brumer 

Chemical Physics Theory Group, 
University of Toronto 
Toronto, Canada M5S 3H6 
(May 30, 2000) 

Extensive coherent control over quantum chaotic diffusion using the kicked rotor model is demon- 
strated and its origin in deviations from random matrix theory is identified. Further, the extent of 
control in the presence of external decoherence is established. The results are relevant to both areas 
of quantum chaos and coherent control. 
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The kicked rotor and its classical limit, the standard 
map, have long served as paradigms for classical and 
quantum chaos The classical dynamics shows char- 
acteristic diffusive energy growth whereas the quantum 
dynamics shows similar chaotic short time behavior, fol- 
lowed by the suppression of diffusion at longer times. In 
this Letter, we demonstrate that the quantum features 
of the chaotic kicked rotor allow for extensive coherent 
control 0] over quantum chaotic diffusion, even in the 
presence of modest decoherence. In particular, we show 
that quantum relaxation dynamics in the kicked rotor 
model is sensitive to the coherence characteristics of the 
initial state, and that altering these characteristics allows 
for control over the energy diffusion. The extent of the 
controlled behavior is vast, from strong suppression to 
strong enhancement of diffusion. 

Consider the kicked rotor whose Hamiltonian is given 

by 



H(L, e,t) = — +X cos(0) S (t/ T - ri ). 



(1) 



where L is the angular momentum operator, 9 is the con- 
jugate angle, I is the moment of inertia, A is the strength 
of the "kicking field" , and T is the time interval between 
kicks. The quantum time evolution operator F for times 
(N - 1/2)T to (N + 1/2)T is given by @ 
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with dimensionless parameters r = HT/I and k — XT/h. 
The classical limit |IJ of this quantum map is given by 
the standard map, which, when expressed in terms of di- 
mensionless variables 9 and the scaled c-number angular 
momentum L = Lt/Tl, takes the following form, 



L 



N 



On — 0N-1 + {Ln + Ln-x)/2 
Ln-i + rcsin(0jv-i + L N -%/2), 



(3) 



where k — kr, and (Ln,9n) represents the phase space 
location of a classical trajectory after N kicks. The sys- 
tem is chaotic for k > n cr = 0.9716.... The resultant dif- 
fusion constant can be defined as the absorption rate of 



the average scaled rotational energy E = (L 2 )/2. Com- 
paring classical and quantum dynamics for typical initial 
classical states shows that quantum dynamics displays 
significant suppression of the classical chaotic diffusion, 
i.e. the external field can only excite a finite number of 
unperturbed energy levels j|. 

The fact that the rotor is a Hamiltonian system and 
the kick is coherent implies that the system maintains its 
quantum phase throughout the evolution. If this is the 
case then the system should be controllable via coherent 
control p|, i.e. by using quantum interference phenom- 
ena to affect the dynamics. To demonstrate this, and to 
examine the extent of possible control, we consider the 
dynamics of states which are initially comprised of su- 
perpositions of two arbitrary angular momentum eigen- 
states, \m) — exp(ir7i0)/v / 27r and \n) — exp(m0)/v / 2 _ 7r. 
Each of these eigenstates is classically allowed, with a 
corresponding classical distribution function given by 



j L/t 



J2-k and p c n {9,L) = 5 L/Tn /2n, respec 



tively g. 

To show that changing the coherent characteristics of 
the initial state significantly alters the subsequent dy- 
namics, we consider the dynamics of states given initially 
by the superposition \ip) = cos(a)|m)+sin(a) exp(i/3)|n). 
Typical results, culled from numerous cases of varying 
a, j3, k and r are shown below and correspond to a weaker 
and stronger chaotic case, and to two values of j3, i.e. 
(3 = 0, and /3 = tt. Specifically, we display below results 
for case (a) \ipt) = (I + 2) ± | - l))/\/2, with r = 0.5, 
k = 5.0 and for case (b) \ipf) = (| + 1) ± | + 2) )/V2, 
with t = 1.0, k = 5.0. Note that neither the basis states 
nor the superposition states are eigenstates of the parity 
operator. 

Figure 1 shows E for each of these two systems and for 
each of the values of (3. Figure la, for example, dis- 
plays E for |^~) ( dashed curve) and for \ip£) (solid 
curve). Clearly, the initial state gives clear diffu- 

sive behavior during the first 40 kicks whereas energy 
absorption in the case of \ipa) ^ s completely suppressed. 
As a result, E(t = 40T) = 9.6 for the \ip~) case, while 
E(t = 40T) = 1.6 for propagation from the initial state 
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\ipa)- Note (a) that this huge difference is achieved solely 
by changing the initial relative phase (3 between the two 
participating states | + 2) and | — 1) in the initial su- 
perposition, and (b) that by contrast, each of | + 2) or 
| — 1) behave similarly to one another with respect to 
energy absorption, giving E(t = 40T) = 5.4 and 6.0, re- 
spectively. Hence, the observed control is due entirely to 
changing the coherent properties of the initial superposi- 
tion state. 
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FIG. 1. The expectation value of the dimensionless scaled 
rotational energy E — (L 2 )t 2 /2h 2 versus time (in units of T). 
(a) Solid and dashed lines are for the initial states \ipt) and 
\tpa), respectively, r = 0.5, k = 5.0. (b) Solid and dashed 
lines are for the initial states \tp^) and \tp^), respectively, 
t = 1.0, k = 5.0. 

Similar control persists for the more chaotic case shown 
in Fig. lb. Here \if>~£~) (solid line) shows extensive chaotic 
diffusion (i.e. compare ordinates scale for Figs, la and 
lb) for up to 45 kicks, giving E(t = 45T) far higher than 
the value of 70.4 and 77.1 reached by propagating either 
of the basis functions | + 1) and | + 2). Further, and by 
contrast, there is essentially no quantum diffusion after 
t = AT for \ip^} (dashed line). Control (not shown) is 
possible for the resonant case as well, e.g., where r = ir/3, 
but it is somewhat less extensive. 

These differences are also reflected in the details of the 
evolving wavefunctions. For example, Fig. 2 shows the 
probability P(m) of finding the system in the state \m) 
at t = 60T. For case (a), P(m) for |m| > 10 is 15.8% and 
3.4% for (3 — 7r and for (3 = 0, respectively. Similarly, 
for case (b) P(m) differs by a factor of 5.4 for the two 
(3 values (3.2% vs. 17.2%) in the probability of exciting 
the rotor to high-energy rotational states |m), |m| > 20. 
In both Figs. 2a and 2b it is evident that the difference 
in final populations resulting from the evolution of the 



two superpositions is an erratic function of m with few 
evident trends. 
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FIG. 2. Probability P(m) of finding the system in the state 
\m) at t = 60T. Results are for the cases shown in Figure 1 

The behavior shown in Fig. 1 is in sharp contrast to 
that which would be observed for the same initial dis- 
tributions propagated classically. These computations 
are shown in Fig. 3 and result from classical propaga- 
tion of the initially non-positive-definite Wigner function 
p w (9,L) associated with the wavefunction cos(a)|m) + 
sin(a) exp(i/3)|n). That is, we classically propagate 

p w (6, L) = cos 2 (a)/4(^ L) + sin 2 (a) p c n (8, L) 

+ — sin(2a) cos(/3 - (to - n)9)S L/T {m+n)/2 . (4) 

for each of \ipf) and \ip^)- In all cases, the classical 
results show [Fig. 3] strong diffusion, characteristic of 
the chaotic dynamics of the standard map. There are 
only small differences in the E diffusion between \ip+) 
and \tpa) an d between \ip£) and |^~). 

Consider then the origins of coherent control of chaotic 
systems in the quantum dynamics, and the behavior 
in the classical limit. To this end we diagonalize the 
quantum map operator F by a unitary operator U, i.e., 



(i\F\j) 



ULU h 



j , where Uij 



(i\U\j),j 



1,2, ■■■ is the eigenvector with cigenphase fa. Af- 
ter the initial superposition state \ip) = cos(a)|m) + 
sin(a) exp(i/3)|n) is kicked N times, we have 
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+ s\n 2 {a)Y J l 2 U* n U*, l U jl U fn e m ^^ 

+[ism(2a)e- < ^Z 2 C7; ro ^ n C/ ii C/;, i e^*-^) 
ijj' 

+c.c] (5) 

where c.c. denotes the complex conjugate of the imme- 
diately preceding term within the brackets. The total 
term in brackets corresponds to interference effects due 
to initial-state coherence. For large N only the j = j' 
terms will survive in the summations due to rapid oscil- 
lations of e lN ^ s ~^i'' . Hence the last two terms reduce 
to l/2sin(2a) e -^E^ 2 ^ \U 3l \ 2 U* m U jn + c.c. 




time (units of T) 

FIG. 3. As in Fig. 1 except that E is calculated by clas- 
sically propagating the initial non-positive-definite Wigner 
function in Eq. (^J). 

If there is no structure in \Uji\ 2 and the eigenvec- 
tor components U* m and Uj n are perfectly indepen- 
dent, as expected from Random-Matrix-Theory then 
53 • \Uji\ 2 U* m Uj n is small and the magnitude of the in- 
terference term is 1 / \f~D times smaller than that of the 
incoherent terms (where D is the effective dimension of 
the Hilbert space) || , interference vanishes and control is 
lost. Hence, the coherent control of quantum chaotic dif- 
fusion relies upon the residual statistical correlations be- 
tween eigenvector components {Ui 3 }. Indeed, in this sys- 
tem the matrix (z|.F|j) is known to display a band struc- 
ture with the bandwidth 2k where the quantity k 2 /N, 
where N is the size of the banded random matrices, pro- 
vides a measure of the statistical deviations from Ran- 
dom Matrix Theory M. Numerical results show that 
k 2 /N < 0.2 is sufficiently small for control to persist. 
This being the case, we obtain a necessary condition to 
ensure the significance of the interference term, namely, 



k < 0.2N k /k or r > k 2 /(0.2N k ), where N K is the mini- 
mum grid size for accurate FFT calculations with r = 1 . 
Numerical studies indicate that, for n < 10.0, N K « 256, 
implying that we require k < 50/k or t > k 2 /50 for 
control. This makes clear that as one approaches the 
classical limit (by increasing k or decreasing r with fixed 
k), coherent control is lost. 

Further evidence that deviations from random matrix 
theory are responsible for control was obtained by exam- 
ining control using a model composed of a banded matrix 
with random matrix elements. Control was obtained in 
this case as well, but was not as extensive as the kicked 
rotor system since the latter displays less random matrix 
character. 

The dynamics of the kicked rotor in the presence of 
decoherence effects has also been examined both exper- 
imentally H and theoretically 0]. The survival of con- 
trol in the presence of decoherence is of interest both 
in general, and for this particular case. To examine 
this issue we introduce a simple decoherence model. 
Here, the quantum map operator between (N — 1/2)T 
and (N + 1/2)T is taken as RF, where R introduces 
random phases into the system. Specifically, R\m) = 
e i27rr£(m,iV)| m ^ where m = 0,±1,±2,---, and £(m,N) 
takes on random values that are distributed uniformly 
between and 1 for each different m or N. Note that 
this model is such that its r — ► 1 limit corresponds to 
measurement-induced quantum diffusion p0| . 

The density matrix p for the dynamics governed by 
RF can then be obtained as an average over many real- 
izations of £(m, N). We take the linear entropy S = Trp 2 
as a useful additional measure of the purity of quantum 
states, and hence of the effect of decoherence. 

Numerical studies show that for r < 0.05, coherent 
control of quantum diffusion is hardly affected by the 
decoherence. For stronger decoherence, e.g. r = 0.15, 
phase control is essentially lost. Examination of the cor- 
responding values oi S(t = GOT) shows that this is consis- 
tent with maintenance of control when the decoherence 
is sufficiently small so that S(t = 60T) > 0.4. Sam- 
ple results are shown in Fig. 4 where we plot the time 
dependence E versus the time in the presence of small de- 
coherence [Fig. 4a, r = 0.05] and of appreciable decoher- 
ence [Fig. 4b, r = 0.15]. Both panels refer to the case of 
stronger chaos [case (b)]. Comparison of Fig. 4a with the 
decoherence-free dynamics in Fig. lb shows that control 
is still significant, but decoherence is beginning to have 
an effect insofar as E at t = 60T is larger in Fig. 4a than 
in Fig. lb. That is, the results show a slight tendency 
towards the classical behavior. In the case of stronger 
decoherence [Fig. 4b] phase control is greatly reduced 
and long time linear diffusive growth of E is observed. A 
careful examination of Fig. 4b suggests that phase con- 
trol persists until t w 20T when the slopes of the dashed 
and solid curve become virtually identical. However, the 
slopes of these curves are still significantly less than those 
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in Fig. 3b, suggesting that quantum coherence is still 
maintained at these longer times. In essence, it appears 
that phase control over the diffusion rate vanishes before 
quantum coherence is completely destroyed. 
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FIG. 4. As in Fig. lb but in the presence of (a) modest 
decoherence and (b) stronger decoherence. 

Note, finally, that the possibility of control does not 
rely heavily on the specific choice of basis states [ [TT|| . 
For example, we have also obtained extensive phase con- 
trol by adding together either negative parity real basis 
states sin(m0) / \pH or by adding together positive parity 
states cosimG) I ypK . Since the states sin(r7i6')/v / 7r satisfy 
the boundary condition of an infinitely deep square-well 
potential V(9) with V(0) = V(2tt) = +oo and since the 
dynamics of a kicked particle in a well is similar to the 
kicked rotor for similar initial states [0], this indicates 
that one can also demonstrate coherent control of chaotic 
diffusion using superpositions of Hamiltonian cigcnstates 
of a kicked particle in a potential well. 

A number of possible experimental demonstrations of 
the proposed control scenario are evident. For example, 
the kicked diatomic molecule Csl |j 13 is a promising 
molecular system for demonstrating controlled quantum 
chaotic diffusion. In this case, preliminary controlled 
laser excitation could be used to prepare the desired ini- 
tial superposition state (which are here superpositions of 
\J,M) and \J',M), where J and J' are the angular mo- 
mentum and M is their projection on the z axis) and 
to vary 0. Alternatively, one can utilize the square well 
analogy described above to experimentally study kicked 
dynamics of a particle in a well. By contrast, implemen- 
tation of atom-optics approach to studies of control 
appears more difficult, insofar as it is necessary to pre- 
pare initial quantum superposition states, a considerable 
extension of previous work Jl5|] . 



In summary, we have shown that quantum chaotic dy- 
namics in systems such as the kicked rotor may be con- 
trolled by varying, in accord with coherent control, the 
coherent characteristics of the initial state. In particular, 
preparing even a simple superposition state comprising 
two basis states allows a wide range of control over the 
diffusive dynamics. The control persists in the presence 
of weak decoherence, but, as expected, is lost with in- 
creasing decoherence strengths. The predicted range of 
control is extensive, providing results of broad theoretical 
and experimental interest. 
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